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Introduction 

Tin*  general  theory  of  chemical  shifts  has  been  well  known  since  the  paper 
of  Ramsey  in  19b0.  With  the  development  of  improved  methods  for  obtaining 
new  Har t ree-Foek  wave  functions  interest  has  again  been  rekindled  to  theoreti- 
cally determine  the  chemical  shift  tensor  in  this  approximation.  Although  the 
final  equations  are  well  established  they  have  derived  by  several  different 
routes  and  often  in  different  formalisms.  Since  one's  ultimate  goal  is  to  not 
only  correctly  predict  the  chemical  shift  but  also  to  understand  the  origins  of 
its  various  contributions,  it  is  clearly  important  to  have  a good  feel  for  how 
the  various  theoretical  terms  arise.  The  purpose  of  the  present  discussion  is 
to  present  what  we  believe  to  be  a thorough  and,  hopefully,  clear  derivation  of 
the  theoretical  expressions  in  what  we  also  believe  to  be  the  most  transparent 
approach.  A variety  of  background  material  is  included  and  we  delve  somewhat 
more  deeply  into  some  of  the  steps  than  would  one  in  a formal  Journal  article. 

Following  a brief  review  of  the  pertinent  Hamiltonian  the  Hellmann-Feymnann 
(2  3) 

theorem  ’ is  used  to  establish  the  general  shift  expression.  This  is  followed 

(4) 

by  a discussion  of  London's  approach  to  the  origin-gauge  problem  and  perturbed 
Hartree-Fock  theory  in  a variable  basis  ;et . Finally,  we  derive  the  chemical 
shift  equations  employing  London's  approach  and  show  how  they  are  trivially 
modified  for  the  case  where  London  orbitals  are  not  used. 

The  Hamiltonian 

In  the  semi-classical  treatment  of  the  interaction  of  a molecular  system 
with  electromagnetic  radiation  the  usual  Born-Oppenheimer  electronic  Hamiltonian 
is  given  by 

- 1 2^  <pi + ! V2  + v ■ vH 


(1) 
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for  a nucleus  N coupled  directly  to  tin  external  magnetic  field  H.  Is  the 

vector  potential  of  the  1**'  electron  while  V contains  all  other  interactions  being 
considered  Including  all  of  the  attractive  and  repulsive  coulomb  interactions. 

The  vector  potential  consists  of  two  parts! 


Af  - HxO  j-R)  + 


V(rrV 


- -K  T 

1 N 


(2) 


The  first  part  is  due  to  the  presence  of  the  external  field  11  and  contains  an 
arbitrary  vector  R due  to  the  fact  that  in  the  selection  of  A^  we  are  free. to 
select  any  gauge  which,  in  this  example,  is  an  arbitrary  reselection  of  the  origin 
of  the  coordinate  system.  No  such  freedom  of  relative  origin  exists  for  the  second 
term,  however,  which  is  the  contribution  to  of  the  nuclear  moment  p^  and  where 
distances  are  measured  with  respect  to  R^,  the  location  of  p^.  The  form 
exhibited  can  further  be  modified  by  the  addition  of  an  arbitrary  term  of  the 

type  grad(f(i  j))  since  curl  (grad(f (r^))»0  and  does  not  change  the  phvsfc.il 
total  magnetic  field.  In  any  complete  calculation  the  choice  of  gauge  is,  like 
the  choice  of  coordinate  system,  immaterial.  Kpsteiu^'*^  has  shown  that  the 
approximate  unrestricted  coupled  Hartree-Koek  theory  is  also  gauge  invariant 
wlthia  any  given  set  of  trial  functions  in  the  sense  of  yielding  identical  energies 
for  an  arbitrary  choice  of  gauge.  The  general  problem  is,  of  course,  that  exact 
calculations  are  not  currently  possible  for  most  molecular  systems  of  interest 
and  liar t ree-Kock  calculations  usually  employ  either  a limited  basis  set  or 
approximations  to  the  necessary  integrals. 

Since  the  direct  nucleus-  field  coupling  is  easily  treated,  the  presence  of  the 
vector  potential  yields  a perturbation  of  the  form  (electron  subscript  and  sum 
repressed) 


F 


^ " 2mc  {A,p+,,‘A)  + 2mc*  A*' 


Sc  A*p  + 2mc?  A*A 


where  the  coulomb  gauge  (V*A"0)  has  boon  selected' 

For  a system  with  multiple  perturbations  such  that 


- J(  4 4 B /('<»>  4 4 ... 


where  a,p  nro  the  usual  ordering  (expansion)  parameters,  wo  may  expand  both  the 
wnvefunotinn  and  energy  in  a similar  series: 


|t|>>  ■ |ij'  > + + ii  1 1J1  ^ ^ > + ap|i/n,P)>  + ... 


c - c + ac(w)  + pt  + aPi^w,(')  + ... 

o 


so  that , for  example | 


(o.p)  ^ ajV. 

3u3P 


c<a>  - It 

do 


a-p -0 


a-p...  . .»q 


The  Mel  1 man- Feynman  theorem 


(2  ,3) 


yields  the  result 


so  that  further  differentiation  leads  to 


a2L  . <<(|32^-|w,s 

3o3P  3p ' 3a  ' ’ v ' Xi  I an  t ’’  l a™ an  I “ 


3a  1 3p 


(10) 
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or,  since  the  mixed  derivative  is  independent  of  the  order  of  differentiation 


32e 
38  da 


< | >+ 

3a  1 38  |V  V|33  '3a  V|3a38|V 


(ID 


a result  which  Moccia  ^ , referring  to  earlier  work  by  Dalgarno  and  Steward^  ^ > 
callsthe  interchange  theorem.  Accordingly,  we  find 

c(a-6)  - <(.<a)l^(B)lv  + l*<a)>  + <*obt<“'e)l*„>  (12) 


where  the  second  term  is  the  hermitian  adjoint  of  the  first  term. 

If  we  write  the  chemical  shift  tensor  a such  that  the  nuclear  spin 
Hamiltonian  is 

^Spln  = - H.(l  - o).pN  . (13) 

then 


3i£ 

“i^j 


(V*vi<)  I.i.(Hj). 


H-v° 


(14) 


and  we  can  calculate  the  general  i-j  element  of  the  tensor  o from  a knowledge 

of  1 1 1>  > and  the  first-order  perturbative  effect  of  the  external  field  H on  |ij<>. 
o 

Use  of  the  interchange  theorem  reverses  the  H and  p symbols  but  must  lead  to  the 
same  physical  answer. 

The  particular  choice  above  of  finding  the  first  order  external-field- 
perturbed  wavefunction  has  the  practical  advantage  that  for  a system  of  many 
nuclei  it  need  be  carried  out  only  once  whereas  had  we  chosen  to  determine 
♦ a calculation  would  be  required  for  each  nucleus  N. 


London  Orbitals 
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The  freedom  to  select  an  arbitrary  gauge  would  seem  to  add  flexibility  to 
the  chemical  shift  calculation  in  terms  of  choosing  a gauge  that  leads  to  a parti- 
cularly transparent  physical  interpretation.  The  fact  that  most  approximate  cal- 
culations are  not  gauge  invariant  diminishes  this  flexibility  and  presents  one 
with  the  problem  of  which  gauge  is  "best"  in  a particular  situation.  The  usually 
trivial  choice  of  coordinate  system  origin  is  now  more  important  since  this  con- 
stitutes one  class  of  gauge  choice.  One  solution  to  this  particular  problem 
is  to  follow  London  and  replace  the  orbital  ^(r)  centered  at  by  what 


we  shall  term  a London  orbital 


(9) 


(LO) 


\,<r) 


^ A(K"),r 

c ♦ (r)  * |n  > 


(15) 


since  now 


-ie 

, . e..  /fie 

(P  + ^jA)Xn  “ e 


A(Rn)*r 


<p  +-(A-A(I^mn 


(16) 


and  the  vector  potential  term  in  the  remaining  operator  becomes 


1 v(r~v 

A-A(Rn)  - ^xCr-R.,.)  4 


(17) 


Now,  in  both  terms  in  the  operator  (applied  to  a particular  orbital  n) , a 

specific  reference  origin  is  referenced:  still  the  nucleus  N for  the  second 

term  and  now  R , the  center  of  orbital  n,  for  the  first  term.  The  use  of 
n 

London  orbitals  thus  suppresses  any  arbitrariness  in  selection  of  origin;  one 
has,  in  fact,  for  each  choice  of  origin  made  a unitary  transformation  of  the 
basis  so  that  so  long  as  one  uses  the  same  basic  orbital  set  ( 4>n } the  gauge 
invariance  is,  as  Epstein ^^puts  it,  "enforced".  Epstein  points  out  that 
this  resulting  gauge  invariance  within  a given  set  of  functions  is  not  import- 
ant for  just  that  point,  but  rather  that  the  physical  significance  lies  in  the 


r 


fact  that  o_lcA(R«0*r^CA  (r)  is  through  first  order  in  H the  exact  eigenfunction 

n 

of  a one-electron,  one  center  uniform  magnetic  field  problem  which  has  as  an 
eigenfunction  when  H=0.  Since  London  orbitals (sometimes  called  GTAO)  account 
for  only  the  translational  invariance,  some  feel  that  the  gauge  problem 

can  be  overcome  only  by  the  use  of  large  basis  sets,  despite  some  recent  succors— 
es^2^obtaincd  by  use  of  the  London  basis. 

Per turbed  Hartree  Fock  Theory  and  Variable  Basis  Sets 

Use  of  London  orbitals  is  an  example  of  working  with  a variable  basis,  that 
is,  a basis  set  in  which  one  or  more  perturbation  parameters  are  explicity  con- 
tained in  the  description  of  the  wave  functions  employed.  Moccia^  has  treated 
this  problem  and  presents  general  formulas  which  are  very  useful  not  only  for 
the  variable  basis  set  problem  but  for  the  more  conventional  approach  in  which 
the  perturbation  ordering  parameters  appear  only  in  the  Hamiltonian.  One  pur- 
pose of  this  section  is  to  briefly  review  the  general  procedure,  to  establish 
notation,  and  finally  to  present  expressions  which  allow  the  determination  of 
presented  before  in  the  general  formula  for  the  chemical  shift  tensor 
(equation  14). 

We  consider  a closed  shell  system  described  by  a set  of  doubly-occupied 
molecular  orbitals  <|^(r)  written  as  a single  Slater  determinant.  The  total 
energy  is  given  by 


? 2fu  + 1 - V 


where  f^  is  the  expectation  value  of  all  the  one-electron  terms  with  ^ and 
Jjj  E (ii/jj)  and  =(ij/ji)  are  the  usual  coulomb  and  exchange  integrals  in 
the  basis: 

(ij  |kl)=/i<<*(r1)i//*(r2)  ~ <l> ^ (r ^ $ t(r 2) dV ^dV 2 


(19) 
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Requiring  E to  be  stationary  with  respect  to  an  arbitrary  change  in  ^ 
subject  to  maintenance  of  the  MO  orthogonality  leads  to  the  Fock  equations 


l <m|f|n>CM  + l l crtCciCn<  <2'mnlrs>  “ (ms | rn) } = l e1jCnj<,l,l«> 


n 


ni  j r,s,n  rj  Sj  ni 


n,j 


(20) 


where  the  MO  has  been  expanded  in  a general  atomic  orbital  basis  |n> 

where 


iv = E ln>c 


ni 


(21) 


and  <m|n>  is  an  overlap  integral  in  the  basis  |n>.  Since  for  this  case 
(doubly  occupied  orbitals) the  density  matrix  P can  be  written 


occ  . 

P ■=  l 2C  .C 
mn  " nj  mj 


(22) 


we  can  rewrite  the  Fock  equation  in  matrix  form  as 


where 


(F  + GtP)}^  = I e^SCj 


mn 


<m | f | n> 


G(P)  = l pRr(<mrUlns>“  T*  mr|g(sn>) 


mn  “ sr 
rs 


S *=  <mln> 
mn  1 


(23) 


(24) 


and  is  the  column  vector  representation  of  in  the  AO  basis  { | n> } . The 
6elf-consistcnt  nature  of  the  approach  is  apparent;  one  assumes  an  initial 
density  matrix  and  iterates  until  self-consistency  is  reached.  Usually,  the 
matrix  e (of  elements  c^)  is  presumed  diagonal.  If  one  is  to  carry  out  a 
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perturbed  Hartree-lock  procedure,  this  may  be  presumed  to  be  so  in  zeroth  order 
but  not  in  higher  orders. 

Following  Moccia  consider  the  first  order  perturbation  equation  in 

for  the  Hartree-Fock  procedure: 


(F°  + G°(P))C.a  + 


(F(a)+  C(a)(P))C* 


, . occ 
o.o.a  , o_(a)_o  c 

e s C + c S C.  + l 


c(a)s°C° 
ij  S j 


(25) 


Here  we  have  used  a diagonalized 
obtain 


ij 


in  zeroth  order. 


By  rearrangement  we 


(F°+G°(P)  - e°S°)C°  + (F(a)+  C(a)(P) 


c°.  S(t°)C* 


?CCe(a)S°C° 
j lj  ‘ J 


(2b) 


Note  that  G(P)  may  contain  a first-order  contribution  from  botli  the  density 
matrix  P and  the  electron  repulsion  terms  <mr|g|ns>  i J_  a variable  basis  set 
is  employed.  By  a procedure  similar  to  ordinary  perturbation  theory,  c”,  may 
be  projected  on  both  the  occupied  and  unoccupied  orbitals  to  obtain 


occ  , . vac 

l < - KV“,CX  + l 


+ cSa\ P) 


(27) 


where,  we  note,  both  the  occupied  and  vacant  (in  zeroth  order)  MO's  enter  into  the 
expression  for  C°,  and  where 


GW(p) 

mn 


l 

r,s 


l P^“^(<mn|g|ns>°-  ~ <mr|g|sn>°) 

+ Psr  ^<mr lslns>(a)  ~ f <mr|g|sn>^a)}] 


(28) 


Since  an  iterative  solution  to  P 


is  desired  whore 


(<0 


occ  , 

p ■ 2 1 cicl 


p(a)=  2 l [ C°(C°)T  + C®(C°)T] 
it  is  useful  to  rewrite  the  expression  for  as 

, OCC  , , N v 


occ 


.o,„orN  + . ~a  ,-Oxt ! 


o of 
vac  C.  C. 
k k 

k (Ei-Ek) 


(F 


c?  s(a>>  C° 

k i 


(29) 


r o _o+ 
\ Ck 


vac 

+ ^ (c°-e°) 

k Vi  V 


G(«)(p)  co 


(30) 


occ  , 

O V O OT 

so  that  in  the  first  term  one  may  recognize  P = 2 £ anc*  *n  terms 

]( 

two  and  three  the  term  C,C  /(e.'-e,)  which  need  be  calculated  only  once,  as 

R K i K 

indeed  is  true  for  the  entire  second  term.  Indeed,  if  we  neglect  the  third 
term  we  have  what  is  referred  to  as  uncoupled  Hartree-Fock  perturbation  theory. 
It  is  the  electrostatic  integrals  that  couple  the  various  C°  equations  that 
require  the  iterative  procedure  in  perturbed  llartree-Fock  and  lead  to  this 
approach  being  referred  to  as  "coupled"  Hartree-Fock. 

While  our  procedure  requires  us  to  determine  the  (or  C^)  in  the  general 

wave-function  expansion 


(13) 


I*!*  El  n >cni  " Hn°>Cni  + ln°>Cni  + ln°>Cni  + "* 
n n 


(31) 


use  of  a variable  basis  introduces  the  perturbation  parameters  in  the  basis  and 
thus  into  a variety  of  integrals.  For  example,  in  employing  London  orbitals 


for  the  chemical  shift  problem 


"^T'A(R.n)*r  =~-  A(R  )»r 

/ (4>me  )*(4>n*c  )dV 


and,  since  A(Rn>=  ^xRn 


S S<P*  4>  ((R  -R  )xr).dV 

mn  J 2fic  ym  yr  m n j 


Similar  terms  enter  in  the  F(“'  and  G(a)  matrices.  The  general  manipulation  of 
such  terms  will  become  more  apparent  as  the  next  section  is  developed. 


The  Chemical  Shift  Expression 

We  wish  to  write  out  explicitely  the  various  contributions  to  the  chemical 
shift  tensor  and  exhibit  their  origin  in  terms  of  the  general  expression  given 
in  equation  14.  We  will  assume  that  London  orbitals  are  to  be  employed  but  will 
be  able  to  see  also  the  resulting  expression  when  ordinary  AO's  are  used 


1.  <V0l#t<"i'"Nj)|V  $ (Il>lj 

The  Simplest  term  to  consider  is  that  in  which  only  the  unperturbed  wave- 
function  enters  and  the  bilinear  H-y  term  arises  in  the  perturbing  Hamiltonian. 


(Hxr) 


**NX  rN 


“2  KH-V'N)(r-rN)  -(Htn)  (PN*r)  ] — y 


where  r^  = r-Rjj.  Accordingly, 


(1,1) 
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If  wo  define  the  vow  vector  Xo  ns  the  gauge lese  basis  vector  whose  components 
are  if  “ ( n>  the  various  rerot  It  order  MO's  k become  X0^  and  1^  becomes 


(ViJ 


2 V Co1  C - Tr  )'  2C°Coi  (X*/?^  ’^X  ) 

1 £ \ Xo^lj  xo  k l k k '‘VWj  V 

" I 


o _i  _/v  (1.1)  |n> 


V <m  -Jl 
urn  1 ‘ L , . 
n,m  1J 


(36) 


We  remind  ourselves  again  that  |n>,|m>  are  gaugeless,  non-1.0  orbitals. 

2.  <^(H1)|  -X(,'Nj)N.o>  + h.n . - (12  + 13)  ^ 

As  we  finally  evaluate  and  1^  to  complete  the  calculation,  the  desirability 
of  dividing  the  Integral  Into  two  parts  will  become  apparent.  Since  t he 
perturbed  wave  function  is  now  involved,  the  perturbed  MO's  x0^'k^*  ^ "*  X "i^C^ 

are  composed  of  two  parts. 


‘b’n  • r «%'y 


Since 


(37) 


occ 

p • I 2ckckT 

k 


P<Hi)  = pd.o) 
i 


ore  (H . ) 

I !«:k<:k  1 


* ‘k<vC  > 


_<•_  (V2j 

<fL'  me  „3 

*N 


(38) 


13. 


we 


finally  obtain  (the  Hcrmitinn  adjoint  now  explicitly  included) 


rMxp) 


(1,).,  - y <P1,0)  <m|— S—  («~ — J — |n> 
2 ij  L nm  i 'me  3 1 


n,nt 
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1^  will  contain  which  contributes  the  part  of  the  H-p  bilinear  term. 

The  derivation  must  be  done  carefully  since  the  perturbation  operator  is  a 
deri  yaUyo  operator  and  will  act  on  the  1^  can  quickly  be  put  in  the  form 
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Since 
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the  contribution  from  the  second  term  takes  on  a different  form  and  we  finally 
obtain 
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where  we  note  the  resemblance  in  form  f the  second  part  of  the  expression  for  1 ^ 

I 

to  that  for  1^  in  equations  36  and  33;  indeed,  they  may  be  combined. 

Expressing  p-  -idiV  we  may  combine  all  terms  to  finally  obtain 
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Having  specified  a particular  perturbation  - the  external  ll-field  and 
internal  p-field-the  perturbation  expressions  may  also  be  cxplicitely  given. 
The  general  matrix  element  for  i'  has  been  given  in  equation  33.  For  1 i 
we  obtain 


<-l«V  V Sr>i  ln> 
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where  t/£co re  is  the  usual  unperturbed  1-election  part  of  the  Hamiltonian. 
The  C^'i^  term  is  given  by  replacing  in  equation  28  P /"'by  P ^'i^where,  for 

16  IS 


example, 
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Equation  A3  represents  the  general  results  for  the  i-j  clement  of  the 
shielding  tensor  calculated  from  the  coupled  perturbed  Hartrce-Fock  approach. 
For  solution  work  we  may  desire  only  the  isotropic  part  of  o obtained  by 
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